Phase directed excitonic transport and its limitations due to environmental influence 
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We investigate theoretically the transfer of excitation along a one dimensional chain of monomers 
for a situation in which initially the excitation is shared coherently by two monomers. We show 
that depending on the relative phase between the two monomers strong directionality of the energy 
transfer is possible. It is also investigated how dephasing, induced by an environment, influences 
this directed transport. 



I. INTRODUCTION 

Controlling the direction of motion of a quantum me- 
chanical wave packed has received much interest, e.g. in 
the context of unidirectional rotational motion* — . In 
the present work we will investigate unidirectional exci- 
tation transfer in assemblies of monomers that interact 
via resonant dipole-dipole interaction (see for instance 
Refs— ~— ). Among the various examples are molecular 
aggregates 1 ^, organic molecular solids^, light harvesting 
systems 1 - , chains of quantum dots, or chains of ultra-cold 
Rydberg atomsi2~— . In all these systems the transport 
of electronic excitation can adequately be described by 
the exciton theory of Frenkel2°r— , where electronic ex- 
citation localized on a monomer can be transferred to 
another monomer due to transition dipole-dipole inter- 
action. 

In most of the studies of this excitation transfer usually 
the initial condition is taken in such a way that the exci- 
tation is localized on a single monomer, i.e. one monomer 
is electronically excited and all the others are in their 
electronic ground state (e.g. in Refs— ' 13 ' 23 ~ — ) Very lit- 
tle attention is paid to the investigation of initial states 
that extend over several monomers. In this situation, 
initial phase relations between different monomers play 
an important role^ 6 -. With the increasing capability to 
create artificial structures, controlling the initial condi- 
tion can be used to tailor the properties of such systems. 
One possibility to create a state where excitation is dis- 
tributed coherently over several monomers, is to use elec- 
tromagnetic radiation that has a definite phase relation 
between different monomers. Since for a typical molecu- 
lar aggregat o 14 ' 16 ! 27 the distance between the monomers 
is in the oder of a few Angstrom and the excitation wave- 
length is in the visible, direct addressing of individual 
chromophores is quite difficult and would probably re- 
quire the use of near-field methods. For chains of ultra- 
cold Rydberg atoms, where the distances between the 
monomers can be in the /an rang o 28 ' 29 , this addressing 
should be easier to achiv and in principle be possible with 
present-day technology. 

In this work it is pointed out that, even for an ini- 
tial state in which the excitation is delocalized only over 
two neighboring monomers, it is possible to obtain very 
strong unidirectional transport along a one-dimensional 
chain by a suitable choice of the relative phase between 



the excitations of neighboring monomers. To investigate 
the robustness of the directionality against fluctuations 
of the surroundings we describe the interaction of the 
Frenkel exciton with its environment within the Haken- 
Reineker-Strobl model, including pure dephasing. 

The paper is organized as follows: In section|TI]we spec- 
ify the system considered and the initial states. In section 
IIIII numerical calculations are performed to illustrate the 
influence of the relative phase on the direction of trans- 
port. Here also the interplay with the environment is 
investigated. In section llVl we give intuitive analytical ex- 
planations for the observed directionality. In particular, 
an exact formula for the time-dependent mean-position 
of the excitation is derived, including the interaction with 
the environment. In section [V] we conclude with a sum- 
mary of our findings and give a brief outlook. 

II. DESCRIPTION OF THE SYSTEM 

We consider a one dimensional array of N monomers. 
The monomers are described by electronic two-level sys- 
tems, where | ) denotes the ground state and | (j)^ ) the 
excited state of monomer n. The energy difference be- 
tween the two states is denoted by e n . Since we focus on 
the transport of a single exctation, only aggregate stats 
of the form 

kn> = |#>---|ft>-H$r>, (!) 

where monomer n is excited and all the others are in their 
ground state, are taken into account. The aggregate is 
described by the Hamiltonian 

N N 

H = ^ e «l 7r ")( 7r « I + X! Vnm\K n ){lZ m \ (2) 
n=l n,m=l 

corresponding to a simple Frenkel exciton mode l 16 ' 30 . 
The transfer matrix element V nm is typically taken to 
be of dipole-dipole form with a distance dependence 
V nm ~ l/\R n — R m \ 3 , where R n and R m are the positions 
of monomer n and m. We focus on a one-dimensional 
regular chain of equidistantly spaced monomers and use 

V nm = V — pr, n ^ m (3) 
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where V scales the interaction strength between the 
monomers. In the following we take e n = e independent 
of the monomer index n. 

Following the approach of Haken, Reineker and 
StroblSiiS, the interaction of the monomers with the sur- 
roundings is taken into account by the time dependent 
Hamiltonian 



H'{t) =2^?n(t)|7r„)<7r B 



(4) 



Here q n (t) are stochastic environment fluctuations which 
couple to the excitation on monomer n. They are taken 
to be real Gaussian random variables with mean 



<<?»(*)> =0 



and variance 



(qn(t)q m (s)) = 7nS nm 6{t 



(5) 



(6) 



The average over the fluctuations is denoted by (• ■ ■ ). In 
Eq. ([6]) we have assumed that fluctuations at different 
monomers n and m are uncorrelated and that the envi- 
ronment has no memory, i.e. it is Markovian. In recent 
years there has been much interest in the influence of 
non-Markovian environments on energy transfer— i^I but 
to discuss the basic influence of the environment on the 
directed transport, a Markovian description will be suffi- 
cient. In the conclusions we will briefly discuss the effects 
of non-Markovian environments together with the influ- 
ence of internal vibrations of the monomers. Further- 
more we will assume that the "f n are site independent, 
i.e. 7„ = 7 for all n. 

From the Liouville-von Neumann equation, averaged 
over the environmental fluctuations, one finds for the ma- 
trix elements p nrn — ( ir n \p\ ir m ) of the averaged density 



matrix 



Pnm \t) 



nm 
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nm) Pnm 
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The probability to find excitation on monomer n is given 
by the population p nn {t). 

In the following, for convenience, we take the number 
of monomers N to be even and assume that initially the 
excitation is located on the two monomers in the middle 
of the chain, i.e. 



PN/2,N/2 
PN/2+1, N/2+1 
PN/2, N/2+1 
PN/2+l,N/2 



PL 
PR 
PLR 

* 

Plr 



(8) 



where p^ denotes the population on monomer N/2 and 
pn that on monomer N/2 + I. The coherence between 
the two monomers is denoted by plr- All other elements 
are zero. In matrix notation this can be written as 



/o 



Pl Plr 
Plr PR 



(9) 



0/ 



We take the initial population on monomer N/2 and 
N/2 + 1 to be equal, i.e. pl = pr = 1/2. Furthermore we 
write for the coherence plr = aexp(— ?0), with a > 
and — 7r < 6 < 7r. Of particular interest will be the phase 
factor between the two monomers. For a pure initial 
state we have a = 1/2 and the corresponding initial wave 
function is given by 



1^ 



A0\ 



(10) 



In the following we show how the phase factor O deter- 
mines the directionality of excitation transfer. 



III. NUMERICAL CALCULATIONS 



In this section we will investigate numerically the time 
development of the populations p„ n (t) by solving Eq. ([7]) 
for different values of the relative phase and the de- 
phasing rate 7. We take the interaction strength |V| as 
unit of energy. Consequently we express time in units of 
h/\V\ and the dephasing rate 7 in units of |V|/ft. 

Exemplarily, in Fig. [JJ the transfer dynamics is shown 
for five values of and two different 7 for a chain of N = 
60 monomers. The coupling V is taken to be positive 
(changing the sign of V leads to a reversal of the direction 
of propagation). Initially the excitation is localized on 
the sites 30 and 31 with a pure initial state as given in 
Eq. (ITU)) . In the left column the fully coherent case 7 = 
is shown. In the right column the case 7 = 0.3 is shown. 
In each panel the time dependence of the mean 



(11) 



is displayed as a black solid line. 

Consider first Fig. [JJc) and h) where = 0. Here one 
clearly sees a symmetric spreading of the initially local- 
ized excitation. The mean stays constant at its initial 
location. This case is similar to the situation where the 
excitation is initially localized only on one site (see e.g. 
Ref^). 

Upon increasing the phase difference to positive val- 
ues the transfer to the left becomes enhanced (while the 
transfer to the right decreases). For a value = +7r/4 
(see Fig. [TJ b,g) already a clear asymmetry can be seen, 
which becomes stronger when increasing 0. For = 
tt/2 the transfer to the left reaches its maximum (see 
Fig.fflai). 
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FIG. 1: Probability p nn (t) to find excitation on monomer n at time t for different values of the phase and the dephasing rate 
7. The black solid line is the time-dependent mean of the distribution p n n{t). Left column 7 = 0, right column 7 = 0.3. Time 
in units of h/V. The values of are indicated in the figures. 



We have here restricted the discussion to the interval 
— 7r/2 < O < 7r/2. The result for an arbitrary phase 
can be simply related to the interval [0,7r/2]. For exam- 
ple, for > 7r/2 we can write = it/2 + 0' and we find 
the same result as for = 7r/2 — 0'. Note that a change 
of sign of just reverses the directions. See Fig. [1] (d,i) 
for = -tt/4 and Fig.[T](e,j) for = -tt/2. 

While for the case 7 = (left column of Fig. []} the ex- 
citation propagates as a well defined "wavefront" over the 
whole time-range of the propagation, in the case 7 = 0.3 
(right column) the exciton dynamics soon becomes simi- 
lar to a diffusion process. This behavior is also reflected 
in the dynamics of the mean. First we note that one 
clearly sees the change of directionality with changing 
phase 0. For the case 7 = the mean moves along the 



chain with constant velocity (characterized by the con- 
stant slope of the black line). This velocity depends on 
0. For a finite dephasing 7 = 0.3 and short times the 
mean also moves with the same velocity as in the case 
7 = 0. Soon, however, it slows down and remains finally 
quasi constant. We will quantify these observations in 
the next section. 

The strong dependence on the phase can be seen 
even better by considering the probability to find excita- 
tion on the left or right side of the chain, i.e. by consid- 
ering the quantities 

N/2 N 

PL® = J2? nn and p «w = J2 p™- ( 12 ) 

«=1 n=£ + l 
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FIG. 2: Asymptotic population on the left side of the chain 
as a function of the initial phase O for the case of vanishing 
dephasing, i.e. 7 = 0. Crosses are the numerical results for 
t = 12h/V , the dashed line is the analytical solution Eq. (|20l) . 
The inset shows the population Pl (t) as a function of time for 
the initial phase <d — tt/2 and different values of the dephasing 
rate 7. From top to bottom the values are 7 = 0.0, 0.1, 0.2, 
0.3 and 0.5. 



In Fig. [2] the probability Pi(i= 12) (i.e. to have popula- 
tion on the left side at the final time t = 12 of Fig. Q}, is 
shown as a function of the initial phase 0. With increas- 
ing initially there is a fast increase in Pl followed by a 
quite fiat region around the optimum at 6 = tt/2 = 90°, 
thus for small deviations AO < 20° from this optimum 
there is still nearly the optimal directionality. The inset 
shows the population Pz,(t) as a function of time for the 
initial phase 9 = tt/2 and different values of the dephas- 
ing rate 7 . 



IV. ANALYTICAL CONSIDERATIONS 

A. The fully coherent case 

In this section we give an intuitive explanation for the 
observed directed transport. Since the important contri- 
bution to directionality stems from the coherent motion 
of the excitation, in the following we will focus on the case 
7 = 0. To obtain simple analytical results we restrict to 
nearest neighbor interaction and impose periodic bound- 
ary conditions (i.e. Vn,i = V). Then the eigenstates of 
Eq. ([2]) are given byi£ 



1 N 



I TTr, 



with 

27T 

k = ^3> j = -N/2,...,N/2-l. 
The corresponding eigenenergies are 
E k = e + 2V cos k. 



(13) 



(14) 



(15) 



With these results we will now discuss how one can under- 
stand directionality from the initial population of eigen- 
states \<f>k)- The sign of the wave vector k appearing in 
the definition of the eigenstates and eigenenergies can be 
viewed as the direction of propagation associated with 
the stationary state Eq. (|13|^ . For N — > 00 the mean 
velocity associated with the wave vector k is given by 
Vk = dkE(k) — — 2Vsin(fc). For instance for k > and 
positive V the velocity is directed to the left. Expanding 
the initial state | ip ml ) defined in Eq. (fTU| w.r.t. the basis 
states Eq. (fT3|) one finds 



1 



E 



-ikN/2 1 



1 



,i(e-k) 



)\4>k). (16) 



The probability to find a certain k within the initial state 
I V> ini ) (defined in Eq. ([ID])) is given by 



1 



p & = K^|V ini )l 2 =-^(i + cos(fc-e)) 



(17) 



From this equation one sees that for = the proba- 
bility to find positive and negative fc-values is equal. For 
negative values of the maximum of the function Pf, is 
shifted to positive fc, i.e. velocities directed to the right 
for V > 0. This has consequences also for the initial 
velocity 



' fe ^ = -Vsin(0), 



(18) 



We 



which shows a sinusoidal dependence on the phase 0. 
will come back to this result in the next section. 

Consider now the probability to find positive fc-values 
in the initial state 



Pfc>o — 



jv/2-: 

E 



~ wJ 2 N 



cot I 



Ni 



N 



In the limit TV — > 00 Eq. (119)) simplifies to 

_ 1 sin(0) 
rk>0 ~ 7T H 

Z 7T 



sin(0). (19) 



(20) 



This analytical function is shown in Fig. [5] as a dashed 
line. It is in very good agreement with the asymptotic 
populations (crosses in Fig. [2]) found numerically for the 
case 7 = 0. Thus Pt>o is just the population on the left 
side of the chain at large times, i.e. P L {t — > 00) = Pfe>o- 
For a finite chain this result holds only for times before 
the first reflection at the end of the chain occurs. Note 
that the time has to be large enough so that the initial 
oscillations shown in the inset of Fig. [5] do not play a role 
anymore. 



B. The first moment of exciton motion 

As seen in Section Mil the mean position of the exci- 
tation M(i), defined in Eq. (fTTj) . is a useful measure to 
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quantify directionality. For our initial condition the ex- 
citation is localized on monomer N/2 and N/2 + 1 with 
equal probability. Thus, M (t) > (N + l)/2 means that 
the excitation is mainly localized on the right side of the 
chain, for M(t) < (N + l)/2 it can mainly be found on 
the left. 

Following the treatment of Reineker— we can find a 
differential equation for M{t). In the following, for sim- 
plicity, we assume nearest neighbor interactions only and 
write V nm = V{5 m ^ n+ i + <5 m .„-i) with equal interaction 
between all neighboring monomers. Furthermore we con- 
sider an infinite chain. Then, after some algebra, one 
finds 

jM{t) = V-m (21) 

with 

4>(t) = i 2j(Pn+l,n - Pn.n+l) (22) 
n 

Differentiating 0(f) w.r.t. time and noting that p nn — 
J2 n Pn+i,n+l and J2 n Pn+2.n = J2 n Pn+l,n-l, from 
Eqs. © and <J22]) we get 

|^)=- 7 0(t). (23) 

This equation can be easily integrated to give 

0(t) = 0(0)^ (24) 

with 0(0) determined by the initial density matrix p mi . 
Note that for our special initial condition Eq. © with 
Plb = a exp(— iO) one has 0(0) = 2 Imp^R = —2a sin O. 
Inserting the result Eq. (|2"4")) into Eq. (j2"Tj) gives 

( V(j){0) -t for 7 = 

M(t) = M(0) + \ (25) 
[ Z0( O )(l- e -7t) for 7 ^0 

From Eq. (j2"5")) one sees that the initial density matrix 
p ml (entering via 0(0)) plays a crucial role in the time 
evolution of M{t). If 70(0) > the mean M(t) will 
move to larger n. Note that a change of the sign of V 
also changes the direction. 

Before considering the initial condition Eq. (j5J) in more 
detail, let us first discuss the dependence on the dephas- 
ing rate 7. From Eq. (|2"5j) one sees that for 7 = 
the mean M(t) moves with a constant velocity given by 
|V0(O)|, which for a pure state (where 0(0) = — sinO) 
is in agreement with Eq. (j!8[) obtained by considering 
the distribution of k- values in the initial state. For 
7^0 this velocity becomes time-dependent. As long 
as t < I/7 one can expand Eq. (|2"5"|) to find M(t) sa 
M(0) + V(f>(0)t + . . . . Thus initially the mean moves with 
the same velocity as without dephasing. For t > 1/7 the 
directionality of the transfer is mostly lost and the mean 
M(t) slowly approaches its long time limit 

M(t > 1/7) = M(0) + -0(0) (26) 



After this discussion of some general properties of M(t) 
we go back to the special initial condition Eq. © with 
Plr = aexp(— iO). Then, as already noted above, one 
has 0(0) = —2a sin 9. For given V and a the direction 
of transport is solely determined by the phase 0. We 
find that for 8 < the transport is directed to the right 
(for positive a and V) with maximum at —n/2 while 
for < < 7r it is directed to the left. For = 
the situation is symmetric and the mean remains M(0). 
Note that a scales the "strength" of directionality For a 
fully mixed state (a = 0) there is no directionality; for a 
pure initial state (a = 1/2) the directionality is maximal. 



V. SUMMARY AND CONCLUSIONS 

In the present work we have investigated how the di- 
rection of excitation transfer along a one-dimensional ag- 
gregate depends on the relative phase between two ini- 
tially excited monomers. It was found that in the case, 
when there is no interaction with the environment, for a 
phase difference of ir/2 more than 80% of the excitation 
propagates into one direction. The direction depends on 
the sign of the phase. For smaller phase differences the 
directionality decreases. The same holds, if the initial 
state is not pure. As expected interaction with an en- 
vironment, which destroys phase relations between ex- 
citations on different monomers, leads to a strong de- 
crease of the directionality. In this case, the direction- 
ality of the propagating wave-packet is only conserved 
on timescales of the order of the dephasing time. This 
dependence on the initial phase is also reflected in the 
time-dependence of the mean of the excitation. For this 
quantity exact analytical results have been derived, show- 
ing in particular that initially the velocity of the mean 
is given by the product of the transfer-interaction be- 
tween the monomers and the sine of the phase-difference 
in the initial state. In the fully coherent case this ve- 
locity remains constant; with dephasing the velocity de- 
creases exponentially on the timescale of the dephasing 
time. While fluctuations of the environment can often be 
helpful to enhance the transport efficiency of excitonic 
energ y 12 i 33 , in the present case the coherence is essen- 
tial for the directed transport to persist. In this context 
it is also appropriate to briefly discuss the influence of 
internal vibrations of the monomers as they are present 
e.g. in aggregates of organic dyes^i. We have performed 
preliminary calculations to obtain a feeling how the di- 
rectionality is affected by these vibrations. For a sin- 
gle undamped vibrational mode per monomer, we used 
the formalism of RefJ^ which is based on the "coherent 
exciton scattering" approximation which allows to treat 
chains up to hundred monomers. We found that the in- 
clusion of a single vibrational mode per monomer barely 
alters the degree of directionality compared to the purely 
electronic case of Section IIV Al However, similar to the 
findings in RefJ^ the velocity of propagation is slowed 
down. To assess the influence of many (damped) vibra- 
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FIG. 3: Engineered interaction between the monomers. 

a) Strength of the interaction between monomer n and n + 1. 

b) Corresponding transfer for an initial state with B = —it/2. 

tional modes we used the non-Markovian Quantum State 
Diffusion (NMQSD) approach of Ref?25., which allows to 
treat continuous spectral densities. For a Markovian en- 
vironment the NMQSD approach can be shown to be 
equivalent to the master equation ([7]). Considering en- 
vironments that are no longer Markovian (e.g. weakly 
damped vibrational modes) the resulting transfer has a 
complicated behavior showing both dephasing as well as 
the influence of individual vibrational modes, which re- 
quires further investigations. However, our investigations 
indicate that the inclusion of vibrational modes does not 
substantially alter the results found in the present article. 

Note that for initial states that consist of coherent su- 
perpositions of a large number of monomers a better di- 
rectionality is possible. For example, we found that when 
the excitation is delocalised over four monomers the di- 
rectionality can already exceed 90%. Clearly, when the 



number of monomers that initially are excited coherently 
becomes even larger a nearly perfect directionality is pos- 
sible, e.g. for a Gaussian distribution of k values around 
a certain mean wave vector kg (this results in a Gaussian 
wave packet in n-space). The present study shows, that 
even for the smallest possible superposition (involving 
excitation of only two monomers) pronounced direction- 
ality is possible. In this regard it is noteworthy that the 
directed transfer also works for arrangements where the 
interaction between neighboring monomers is no longer 
equal. In nature this is the case e.g. in the LH2 com- 
plex of purple bacteria. In artificial systems one might 
design the couplings (e.g. by changing the distance or 
orientations between the monomers) such that the prop- 
agating wave packet has certain properties like a minimal 
spread over the monomers or the refocusing at a certain 
monomer. That such tailoring of the aggregate proper- 
ties can be combined with the directionality induced by 
the initial phase is exemplarily shown in Fig. |3]for near- 
est neighbor interactions as given in panel (a)2I, cho- 
sen such that the initial excitation is focused at the end 
points of the chain. This might stimulate new ideas how 
to employ the excitonic phase in molecular devices. The 
directionality imposed by the relative phase in coherent 
superpositions of localized states might even play a role in 
the energy transfer in natural light harvesting systems 35 , 
when partly delocalized excitations "hop" from one com- 
plex to the next and the exciton remains delocalized over 
a small number of monomers. 
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